LETTER 


doi:10.1038/naturel3304 
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The explosive fate of massive Wolf-Rayet stars 1 (WRSs) is a key 
open question in stellar physics. An appealing option is that hydro- 
gen-deficient WRSs are the progenitors of some hydrogen-poor 
supernova explosions of types lib, lb and Ic (ref. 2). A blue object, 
having luminosity and colours consistent with those of some WRSs, 
has recently been identified in pre-explosion images at the location 
of a supernova of type lb (ref. 3), but has not yet been conclusively 
determined to have been the progenitor. Similar work has so far 
only resulted in non-detections 4 . Comparison of early photometric 
observations of type Ic supernovae with theoretical models suggests 
that the progenitor stars had radii of less than 10 12 centimetres, as 
expected for some WRSs 5 . The signature of WRSs, their emission 
line spectra, cannot be probed by such studies. Here we report the 
detection of strong emission lines in a spectrum of type lib super- 
nova 2013cu (iPTF13ast) obtained approximately 15.5 hours after 
explosion (by ‘flash spectroscopy 1 , which captures the effects of the 
supernova explosion shock breakout flash on material surrounding 
the progenitor star). We identify Wolf-Rayet-like wind signatures, 
suggesting a progenitor of the WN(h) subclass (those WRSs with 
winds dominated by helium and nitrogen, with traces of hydrogen). 
The extent of this dense wind may indicate increased mass loss from 
the progenitor shortly before its explosion, consistent with recent 
theoretical predictions 6 . 

Wolf-Rayet stars are massive stars stripped of their outer, hydrogen- 
rich envelope. These stars blow strong, hydrogen-poor winds. The inner 
part of the wind engulfing the star is dense and optically thick, and 
efficiently absorbs the ionizing continuum from the hot stellar surface. 
Farther from the star, the density drops and the wind becomes optically 
thin in the continuum, leading to a rich emission spectrum of recom- 
bination lines. Detailed models of such spectra can be constructed 1,7 and 
depend essentially on only three parameters: the effective temperature, 
T e ff, in the line-forming region, a normalized radius, R t (a combination 
of the stellar radius, luminosity, mass loss rate and wind terminal 
velocity 7 ), and the chemical composition, Z, of the wind (assumed to 
be uniform, spherical and of a constant mass loss rate). The composition 
of the wind determines Wolf-Rayet spectral classes: stars with dominant 
He and N lines belong to the WN class (with those also showing traces 
of H usually denoted as WN(h)), stars with strong carbon lines belong 
to the WC class and rare (and possibly hotter) stars with oxygen-rich 
spectra belong to the WO class 1 . 

Shortly after a WRS explodes as a supernova, the outer parts of the 
wind (which in some cases 1 extend to radii of > 10 13 cm) that have not 
yet been swept up by the expanding supernova ejecta will emit strong 
recombination lines in response to ionizing flux released by the explosion 
shock breakout from the stellar surface. We estimate an increase in ion- 
izing luminosity by a factor of order 1 0 2 - 1 0 4 , assuming an initial absolute 
magnitude range of —2.5 mag < M < — 10 mag for the exploding WRS 1 


and a typical early- time luminosity of M = — 12.5 mag for the resulting 
supernova 3,5,8 . The effective temperature, T e g, will also change, being 
very high (>10 5 K) shortly after explosion 9 and decreasing as the 
shocked supernova ejecta cool. The radiation illuminating the surviv- 
ing Wolf-Rayet wind will thus effectively vary monotonically through 
the range of temperatures in WRS line-forming regions. Because the 
wind parameters (composition, mass loss rate and terminal velocity) 
do not change, the measured line spectrum observed shortly after explo- 
sion should be similar to that of a WRS with the spectral class of the 
exploding progenitor (because the spectral classes mainly reflect the 
wind composition). The high wind densities around WRSs (with electron 
densities of n e = 10 n -10 12 cm 3 ) imply short recombination times 10 , 
f r ec~ 3.9 X 10 12 (n e /cm 3 ) '(T/K/K) 0 ' 85 s, typically of a few minutes 
for Wolf-Rayet densities and temperatures, and so the emitted spectrum 
will promptly react to the rapidly evolving supernova radiation field. 

We obtained rapid spectroscopic observations of the recent type lib 
supernova SN 2013cu (iPTF13ast) shortly after shock breakout (by 
flash spectroscopy; see Methods). This event was first detected by 
the Intermediate Palomar Transient Factory (iPTF) survey 11 on 2013 
May 3.18 (we express dates in coordinated universal time (utc)), photo- 
metrically confirmed 5.8 h later and promptly identified by an on-duty 
astronomer who triggered rapid follow-up observations 12 , including an 
optical spectrum obtained just 4 h later. Analysis of the early- time light 
curve of this supernova (Extended Data Fig. 1) suggests that it exploded 
on May 2.93, implying that the first iPTF detection and the first spec- 
trum correspond to only 5.7 h and 15.5 h after the explosion, respect- 
ively. A full description of the supernova and its evolution will be 
reported in a forthcoming publication (A.G.-Y. et al, manuscript in 
preparation). We note that this event was independently observed by 
the MASTER survey on May 5.3 (—2.3 d after explosion), and it was 
assigned the name SN 201 3cu following spectroscopic confirmation 13 . 

Our first spectrum of SN 2013cu reveals a continuum and emission 
lines that bear a striking resemblance to spectra of WRSs (Fig. la). 
According to accepted Wolf-Rayet terminology 1 , the spectrum is clas- 
sified as subclass WN6(h) (Fig. la, blue trace); the relative strength of 
nitrogen to carbon lines precludes a WC classification, and the absence 
of any high-excitation oxygen lines is inconsistent with a WO star. The 
stronger lines (Ha, Hp, Niv 17,115 (7, 115 A wavelength) and Hell 
15,4 1 1 ) exhibit a complex profile (Fig. 2) consisting of a relatively broad 
base (~ 2,500 km s 1 full width at zero intensity (FWZI)) on which 
prominent narrow, unresolved lines (FWZI — 3 A; velocity dispersion, 
<150 km s : ) are superimposed. This is consistent with predictions 
for Wolf-Rayet pre-supernova wind velocities 14,15 , although we cannot 
exclude the possibility that at least some of the observed line broadening 
is produced by electron scattering rather than genuine velocity disper- 
sion. To the best of our knowledge, no similar spectra of a stripped 
(H-poor) supernova have been acquired previously. Wolf-Rayet-like 
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spectroscopic features have been observed in spectra of some H-rich 
(non-stripped) supernovae obtained at substantially later epochs, and 
their typically much broader lines were interpreted as emerging from 
interaction with circumstellar material 16 (CSM). We further discuss these 
previous observations in Extended Data Fig. 3. 

We analyse our very early spectrum using the PoWR grid of Wolf- 
Rayet spectral models 7 (http://www.astro.physik.uni-potsdam.de/~wrh/ 
PoWR/powrgridl.html). We find an excellent fit with WN(h) models 
calculated assuming an H fraction of 20% (by mass) and a temperature of 


Figure 1 | Spectroscopy of the type lib supernova SN 2013cu reveals 
transient Wolf-Rayet-like features, a, The early spectrum of SN 2013cu 
(black; specific flux, F versus wavelength) is compared with models of 
WNL-class WRSs 7 (red and grey), showing remarkable similarity in continuum 
shape and in line features (strong He, N and Balmer lines indicate a WN6(h) 
classification; compare with HD 192163 (blue)), b, Emission line evolution 
during the first week. The first spectrum (red) is compared with later spectra. 
By day 3 (blue and cyan), the initially strong Wolf- Rayet features disappear, but 
the Ha line remains constant at 3.4 X 10 -15 ergs -1 cm -2 . The spectrum on 
day 6 (magenta) is almost featureless, except for weak Ha emission (inset) with 
an intensity less than a tenth that of day 3, probably because the line-forming 
region has been cleared by the expanding ejecta, c, SN 2013cu is of type lib. 

A spectrum of SN 2013cu 69 days after explosion (black) is compared with 
the prototypical typellb supernovae SN 1993J (60 days; red) and SN 201 ldh 
(43 days; green), and with the typical non-stripped typell-P supernova 
SN 2004et (45 days; blue). SN 2013cu, SN 2011dh and SN 1993J exhibit strong 
blue-shifted Hei absorption at 5,876, 6,678, 7,065 and 7,281 A (marked with 
black vertical lines), which are not detected in type II-P supernovae. See 
Methods for more details. 

around 50,000 K. This temperature is consistent with the lower limit 
obtained from early ultraviolet photometry from NASA’s Swift space- 
craft (Extended Data Fig. 1). The essentially perfect match of the ob- 
served and modelled continuum shapes indicates that dust reddening 
must be negligible; any pre-existing dust must have been destroyed by 
the supernova explosion flash (Methods). We note that among the large 
catalogue of Galactic WN stars specifically modelled in this manner 7 , 
no stars drive winds that require more than 56% H by mass. Presum- 
ably, custom spectral fits 7 could be calculated and used to determine 
more accurately the physical parameters of the detected Wolf- Rayet 
wind. Assuming that the spectrum was obtained 15.5 h after explosion 
and that a standard ejecta velocity of 1 0 4 km s 1 is applicable, the narrow- 
line-emitting material must be located at radii greater than ~ 5 X 10 13 cm 
for it not to have been swept up and accelerated by the expanding ejecta. 
This lower limit is consistent with the extent of some Wolf-Rayet 
winds, where the line- formation region extends out to several (five to 
ten) times 17 the hydrostatic radius. Recent pre-supernova Wolf-Rayet 
models 14 suggest hydrostatic radii of 10-20 solar radii for WN super- 
nova progenitors, consistent with line-formation regions extending to 
several hundred solar radii 1 , or >10 13 cm. 

We further constrain the physical location of the wind using the fol- 
lowing calculation. We measure the Ha line flux from our spectrum 
(calibrated to our host-subtracted photometry) and obtain f Ho( = 3.4 X 
10 15 ergs 'em 2 . We translate this to line luminosity using a lumin- 
osity distance to the host galaxy, UGC 9379, of d = 108 Mpc (calculated 
for a flat A cold dark matter cosmology with Hubble constant H 0 = 
73 km s 1 Mpc ', a matter density of Q m = 0.27 and a redshift z = 
0.025734, obtained from the NASA Extragalactic Database (NED), as 
well as negligible extinction), obtaining L H[X = 4.8 X 10 39 ergs ’. We 
can then estimate 18 the pre-explosion H mass loss rate 

M = 0.01 (i H a/(2 x 10 39 erg s -1 )) 1/2 

x (v w /(500 km s -1 )) (r/(10 15 cm))^ 2 M 0 yr -1 

assuming that the lines are formed at a radius r by recombination, that 
the width of the emitting shell is similar to its radius, that the explosion 
is spherical symmetric and that the density of the wind decreases as r 2 . 
We assume a wind velocity of v w = 2,500 km s _1 , consistent with our 
spectra and as expected for WRSs, but our results are not sensitive to 
this value and remain essentially unchanged for velocities in the range 
100 < v w < 2,500 km s '. 

We check for self-consistency by calculating the electron density 
and, hence, the Thomson optical depth 

T = 0.3(M/(0.01M Q yr -1 )) 

x (v w /(500 km s -1 )) 1 (r/(10 15 cm)) 1 
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Figure 2 | Emission line velocity structure at 15.5 h. The strongest lines extending out to ~2,500 km s . Other weaker lines are narrow 

(Ha, Hp, Hen 25,411 and the Niv 27,115 complex) show broad wings and unresolved. 


expected at this radius, and require it to be lower than t = 1 for the line 
emission to escape. W e find that this self-consistency requirement places 
a lower limit of r = 2 X 10 14 cm on the radius of the line-formation region, 
and implies substantial mass loss rates, Mw 0.03 Mg yr _1 . If we inter- 
pret the disappearance of essentially all emission lines from our day-6 
spectrum (Fig. lb) as evidence that the wind was swept up by the expand- 
ing ejecta (moving at 10 4 kms *), the radius of the line-emitting region 
must be r < 5.2 X 10 14 cm, which is fully consistent with our estimates. 

We can then calculate the total H mass by integrating over r: 

M tot = 0.006 (m/(0.01M o yr _1 )) 

x (v w /(500 km s -1 )) 1 (r/(10 15 cm)) 1 Mq 

This indicates a range of 0.008MqV M tot v 0.0035Mq for the range of 
permitted H masses. Assuming that the typical H abundance for WN (h) 
stars (—20%) applies, the total wind mass (dominated by He) can be 
estimated to be several times larger than these values. Detailed simula- 
tions 19 show that as little as 0. IMq of He-dominated CSM would result 
in strong spectroscopic interaction signatures (that we do not observe), 
consistent with our derived total masses. 

We conclude that we have directiy detected a Wolf-Rayet-like wind 
from the supernova progenitor with a WN(h) spectral class, indicating 
a low H mass fraction. Assuming that the wind composition we measure 
represents the surface composition of the progenitor star, our observa- 
tions indicate that some members of the spectroscopic WN(h) Wolf- 
Rayet class explode after having lost most of their hydrogen envelope, 
exposing the CNO-processed, N-rich He layer below. Analysis of 
photometric and spectroscopic follow-up observations (A.G.-Y. et al, 
manuscript in preparation) confirms that the explosion was indeed a 


supernova of typellb (Fig. lc), as expected if the progenitor was a 
massive star that lost all but ~ 0. IMq of its H envelope 20 . 

Our observations have interesting implications. First, we note that 
the derived values of the mass loss rate and emission-line-region size 
are quite extreme compared with known Wolf-Rayet observations and 
radiatively driven models 21 , including models with clumpy, inflated 
atmospheres 22 . This suggests that the mass loss rate from the progen- 
itor star may have increased shortly (of order 1 yr for the assumed velo- 
cities) before its explosion. Interestingly, such pre-supernova activity may 
be explained by recent wave-driven models 6 , or a more extreme envelope 
inflation 22 may be indicated. These data can thus provide a key diagnostic 
of the final stages of nuclear core burning in massive stars, which are 
currently poorly understood, with possible implications for the explo- 
sion mechanism itself. In any case, the star probably exploded inside a 
thick wind, and the explosion shock may have broken out from the 
opaque inner wind rather than from the hydrostatic surface of the star 9 . 

Our finding is in general accord with some previous work on type lib 
supernova progenitors. Direct imaging of the progenitor of SN 2008ax 
(ref. 23) is consistent with a WN(h) progenitor. Furthermore, increased 
mass loss during the final year before explosion may inflate the appar- 
ent photospheric radius of the pre-supernova star, making stars with 
compact cores appear to have extended (low-mass) envelopes 22 and pos- 
sibly reconciling the conflicting findings about the progenitor of the 
type lib supernova SN 201 ldh (refs 24-27). Regardless of the exact mech- 
anism, our observations suggest that substantial Wolf-Rayet-like winds 
pre-date at least some type lib supernovae. A strong metallicity depen- 
dence of this process may explain the trend in the type Ilb/type lb 
supernova number ratio with host-galaxy metallicity 28 . Future studies 
of numerous additional supernova progenitors via their spectroscopic 
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wind signatures (Methods) would provide powerful constraints on the 
final stages of massive-star evolution. 

METHODS SUMMARY 

Photometry. The iPTF survey telescope was used to obtain r band observations. 
Photometry was measured using our custom pipeline performing point spread func- 
tion (PSF) photometry on iPTF images after removing a reference image con- 
structed from pre-explosion data using image subtraction. Swift ultraviolet absolute 
AB magnitudes (Extended Data Fig. 1 ) were measured using standard pipeline reduc- 
tion and were corrected for host-galaxy contamination using late-time Swift images. 
Spectroscopy. Our earliest (15.5 h) and latest (69 d) spectra were obtained using 
the DEIMOS spectrograph mounted on the Keckll 10-m telescope using the grat- 
ing with 600 lines per millimetre and an exposure time of 600 s. Additional spectra 
were obtained using ALFOSC mounted on the 2.56-m NOT telescope, LRS mounted 
on the 10.4-m HET telescope (day 3), and LRIS mounted on the KeckI 10-m tele- 
scope (day 6). All spectra were reduced using standard pipelines and are digitally 
available on WISeREP 29 . The method of flash spectroscopy is described in detail 
in Methods. 

Online Content Any additional Methods, Extended Data display items and Source 
Data are available in the online version of the paper; references unique to these 
sections appear only in the online paper. 
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Photometry, r-band images were obtained by the iPTF survey camera mounted 
on the Palomar 48 -inch Schmidt telescope 11,30 . Photometry was measured using 
our custom pipeline performing PSF photometry on iPTF images after removing a 
reference image constructed from pre-explosion data using image subtraction. 
Swiff ultraviolet absolute AB magnitudes (Extended Data Fig. 1) were measured 
using standard pipeline reduction 31 and were corrected for host-galaxy contam- 
ination using late-time Swift images. 

Spectroscopy. Our earliest (15.5-h) and latest (69-d) spectra were obtained using 
the Deep Imaging Multi-Object Spectrograph (DEIMOS 32 ) spectrograph mounted 
on the Keck II 10-m telescope using the 600 line mm 1 grating and an exposure 
time of 600 s. Additional spectra were obtained using ALFOSC mounted on the 
2.56-m NOT telescope, LRS mounted on the 10.4-m HET telescope (day 3) and the 
Low Resolution Imaging Spectrometer 33 (LRIS) mounted on the Keck I 10-m 
telescope (day 6). All spectra were reduced using standard pipelines and are digit- 
ally available on WISeREP 29 . 

Spectroscopic observations are presented in Fig. 1. In Fig. la, the early spectrum 
of SN2013cu reveals Wolf-Rayet wind features. The spectrum (black) is com- 
pared with WNL models 7 (red and grey curves, offset vertically for clarity) showing 
remarkable similarity, both in line features (major species marked; strong He and 
N lines accompanied by Balmer lines indicate a WN6(h) classification) and in the 
continuum shape (demonstrated by overplotting the 56-kK model on the spectrum). 
The similarity in continuum shape to hot model spectra constrains any dust red- 
dening to be minimal, indicating that any pre-existing circumstellar dust must 
have been destroyed; compare with the observed spectrum of the WN6(h) star 
HD 192163 (blue). Consistent with this conclusion, we detect no trace of Na D 
absorption lines. Figure lb shows the emission line evolution during the first week. 
We replot the first spectrum (red), compared with spectra obtained ~3 d after explo- 
sion (using ALFOSC (blue) and LRS (cyan)) and ~6 d after explosion (magenta; 
additional details in A.G.-Y. et al, manuscript in preparation). Spectra were scaled 
and offset for clarity with respect to the high-quality Keck/LRIS spectrum. Con- 
tinuum shapes are identical (within calibration uncertainties) and consistent with 
the Rayleigh-Jeans tail of a hot Planck curve, indicating that the black-body peak 
remains in the ultraviolet at 6 d past explosion. By day 3, the initially strong Wolf- 
Rayet features disappear, with the exception of weaker He n fines and the C m/N hi 
complex around 4,640 A (marked), whereas the Ha line remains constant at 3.4 X 
10 15 ergs -1 cm -2 (with decreased equivalent width due to higher continuum level). 
The spectrum with a high signal-to-noise ratio on day 6 is almost featureless, except 
for weak Ha emission (inset) with an intensity less than a tenth that of day 3, pro- 
bably because the line-forming region has been cleared by the expanding ejecta. 
Figure lc shows evidence that SN 2013cu is of typellb. A spectrum of SN 2013cu 
69 d after explosion (black) is compared with the prototypical type lib supernovae 
SN 1993J at age 60 d (red) and SN 201 ldh at age 43 d (green), and with the typical 
non- stripped type II -P supernova SN 2004et at age 45 d (blue). To allow for slight age 
differences and expansion-velocity variations, we align all spectra in wavelength 
using the weak blue-shifted Fe ii lines at 4,924, 5,018 and 5, 169 A (marked with black 
vertical fines), because they are good tracers of the photosphere. Supernovae SN 
2013cu, SN 201 ldh and SN 1993J exhibit strong He I absorption at 5,876, 7,065 and 
7,281 A (marked with black vertical lines), whereas the weaker 6,678 A absorption is 
similar in SN 2013cu and SN 201 ldh. These He i lines are not detected in the type II- 
P supernova spectrum at all. However, SN 2013cu shows weaker Balmer absorption 
(marked with magenta vertical fines), and the Ha absorption is not clearly defined. 
Both the spectral similarity to supemovae SN 1993J and SN 201 ldh and the strong He i 
fines relative to H indicate that SN 2013cu is spectroscopically a type lib supernova. 
Flash spectroscopy. We define as ‘flash spectroscopy’ the technique of obtaining 
a set of spectroscopic data shortly enough after a supernova explosion that the 
observed spectrum is dominated by features related directly to the effects of the 


shock breakout flash. In particular, flash-ionized CSM recombines and forms 
strong emission fines, revealing, for example, the elemental abundance and, thus, 
the Wolf-Rayet class of a supernova progenitor. In addition, emission fine spectra 
provide a handle on the early temperature evolution, which is difficult to measure 
even using Swift ultraviolet photometry because the black-body peak is initially too 
far into the ultraviolet. This study provides strong motivation for future investi- 
gations using dedicated rapid-response spectrographs, such as FLOYDS 3,34 and 
SEDM 35 , responding to real-time triggers from high-cadence wide-field surveys 12 . 

Although Wolf-Rayet supernova progenitor stars are difficult to study using 
pre-explosion imaging (owing both to intrinsic low luminosity in the optical and 
infrared bands 14,36 and to the possible confusing effect of a bright O/B companion 3 ), 
we demonstrate that they are amenable to study using the flash spectroscopy 
method. WRSs belonging to the WN(h) class may have the most extensive 
winds 1,14 . Application of this method to WC and WO stars may require flash spec- 
troscopy at even earlier epochs (~ 1 h after explosion), before the supernova ejecta 
sweep up the high-density wind. Reducing the latency between the explosion 
and the spectroscopy by an order of magnitude relative to our observations of 
SN2013cu is possible using recently commissioned instrumentation (Extended 
Data Fig. 2). 

Unlike studies of supernova progenitors through pre-explosion imaging, the 
flash spectroscopy method can be applied to relatively distant objects (SN 2013cu 
is located 108 Mpc away, well beyond the 20-Mpc distance typical for pre-explosion 
studies) and to events in galaxies having no high-quality pre-explosion imaging, 
such as a large population of little-studied dwarf galaxies. Judging from local super- 
nova rate measurements 37 , ~300 explosions occur within 100 Mpc of Earth every 
year and can be potentially studied in this manner. The method thus allows routine 
spectroscopic studies of supernova progenitors that were previously only possible 
by extreme serendipity (for example for the progenitor of SN 1987A; ref. 38). 
Within a few years, there will be enough data for the flash spectroscopy method 
to be used to chart wind signatures from numerous supernova progenitors; in 
particular, it will be possible to study systematically the Wolf-Rayet progenitor 
population of stripped supernovae. We thus expect that this method will be 
broadly applied in the coming years. 
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Extended Data Figure 1 | The r-band light curve of SN 2013cu. A parabolic 
model of the flux-time (red solid curve) describes the pre-peak data (lcr 
error bars) very well. Backward extrapolation indicates an explosion date of 
UTC2013 May 2.93 ± 0.11 (MJD = 56414.93; 5.7h before the first iPTF 
detection; see inset); we estimate the uncertainty from the scatter generated by 


modifying the subset of points used in the fit. Our first Keck spectrum was 
obtained about 15.5 h after explosion (vertical dotted line). Early Swift 
ultraviolet photometry (diamonds) places a lower limit of T = 25,000 K on the 
black-body temperature measured 40 h after explosion. 
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Extended Data Figure 2 Flash spectroscopy: rapid spectroscopic 
observations of supernovae during or shortly after shock breakout. This 
simulated SEDM spectrum (red) created by downgrading the observed Keck 
spectrum (black; resolution, R = 2,000) to the coarse SEDM resolution 
(. R = 100) shows that the strong Wolf-Rayet lines (in this case the marked H, 
He and N lines) are still easily detectable and allow us to determine the Wolf- 
Rayet spectroscopic class. The SEDM is an IFU low-resolution spectrograph 
designed for robotic response to transient events, to be mounted on the 
Palomar 60-inch telescope almost continually. Responding to real-time triggers 


from the iPTF survey operating on the same mountain, this instrument should 
be able to obtain low-resolution spectra within ~ 1 h of object detection. 
Because SEDM operates on a smaller telescope than Keck, SEDM data of 
similar quality to the simulated spectrum will require a relatively long 
integration. However, SEDM will be able to observe objects with much reduced 
latency, thus benefiting from stronger line intensities expected from the 
stronger shock breakout flash luminosity processed by a denser wind close to 
the progenitor star, potentially compensating for its reduced absolute 
sensitivity. 
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Extended Data Figure 3 Comparison with early ‘Wolf-Rayet’ spectra of 
SN 1998S. Wolf-Rayet-like features similar to those we observed were 
previously noted in two cases, SN 1983K (refs 39,40) and SN 1998S (ref. 16), 
and persisted for many days after explosion. The spectra of SN 1 983K, classified 
as a type II-P supernova, are unfortunately not available for comparison. 
Spectra of SN 1998S (type Iln) are shown here. The spectra have a similar blue 
continuum slope and a similar Ha profile. The Hell ,1.5,411 and the N IV 1.7, 1 1 5 
complex are weaker in SN 1998S, and the strong lines of N and C (5,806 A) 
are broad, consistent with an origin in shocked CSM 16 rather than in an 


undisturbed Wolf-Rayet wind. The inset shows a close-up view of the strong 
N iv 77,115 complex. N emission from SN 1998S is weak compared with 
He I 77,065 and shows a smooth, broad profile, whereas SN 2013cu exhibits 
a broad base (full- width at half-maximum, 2,500 km s ’) as well as strong and 
narrow (unresolved) N IV lines. These observations are consistent with an 
origin for the N emission of SN 1998S in shocked (and perhaps N-rich) CSM 16 , 
whereas the narrow He I lines may come from a more distant, photo-ionized 
wind. In SN 2013cu even the narrow N lines are much stronger than He, 
indicating a wind with a Wolf-Rayet-like composition at all radii. 
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